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Abstract. In the setting of an n-dimensional Euclidean space, the du- 
ality between velocity fields on the class of admissible bodies and Cauchy 
fluxes is studied using tools from geometric measure theory. A gener- 
alized Cauchy flux theory is obtained for sets whose measure theoretic 
boundaries may be as irregular as flat (n — I)-chains. Initially, bodies 
are modeled as normal n-currents induced by sets of finite perimeter. 
A configuration space comprising Lipschitz embeddings induces virtual 
velocities given by locally Lipschitz mappings. A Cauchy flux is defined 
as a real valued function on the Cartesian product of in — I )-currents 
and locally Lipschitz mappings. A version of Cauchy's postulates im- 
plies that a Cauchy flux may be uniquely extended to an n-tuple of flat 
(n — I )-cochains. Thus, the class of admissible bodies is extended to 
include flat n-chains and a generalized form of the principle of virtual 
power is presented. Wolfe's representation theorem for flat cochains 
enables the identification of stress as an n-tuple of flat (n — 1 )-forms. 



1. Introduction 

This work presents a framework for the formulation of some fundamental 
notions of continuum mechanics. Specifically, using elements from geometric 
measure and integration theory, we consider, within the geometric setting 
of R n , the class of admissible bodies, configurations of bodies in space, the 
configuration space, virtual velocities, and stress theory. 

Cauchy's stress theorem is one of the central results in continuum mechan- 
ics. It asserts the existence of the stress tensor which determines the traction 
fields on the boundaries of the various bodies. As the traditional proof relies 
on locality and regularity assumptions, from both the validity and the appli- 
cability aspects, stress theory is closely associated with the proper choice of 
the class of bodies. Furthermore, an appropriate class of bodies should allow 
the formulation of the Gauss-Green theorem or a generalization thereof. 

In light of these observations, formulations of the fundamentals of contin- 
uum mechanics have considered, since the middle of the 20th century, the 
appropriate choice of the class of bodies. In [22], Noll defined a body as 
a differentiable three-dimensional compact manifold with piecewise smooth 
boundary that can be covered by a single chart and is endowed with a mea- 
sure space structure. The configurations of the body in space provide charts 
on the body manifold and a part of the body is defined as a compact sub- 
set of the body with piecewise smooth boundary. In [38| p. 466] Truesdell 
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and Toupin ignore the formal issue of admissible bodies and tacitly assume 
smoothness wherever necessary. Later on, in [37} p. 4], Truesdell adopted 
the structure suggested in Noll's [22J. The common ground for these early 
works is in the assumption that bodies in continuum mechanics should have 
a smooth structure so that the classical versions of the notions of mathemat- 
ical analysis apply. In |22| . it is shown that Cauchy's original postulate on 
the dependence of the traction on the exterior normal may be replaced by 
an additivity assumption on the system of forces and the principle of linear 
momentum. 

Gurtin and Martins introduce in [10J the notion of a Cauchy flux in order 
to represent the collection of total forces applied to the collection of surface 
elements. A Cauchy flux is defined as an additive, area bounded set function 
acting on the collection of compatible surface elements of the body, and a 
weakly balanced Cauchy flux is defined as a volume bounded Cauchy flux. 

It seems that pQ and |43] were the first to propose that the class admissible 
bodies in continuum physics should consist of sets of finite perimeter. In 
Ziemer's work, admissible bodies are defined as sets of finite perimeter and 
a weakly balanced Cauchy flux is shown to be represented by a measurable 
vector field. The works I24| . which followed, further extended these 
studies. In the class of admissible bodies is defined as the class of 

normalized sets of finite perimeter while in |24j . admissible bodies are defined 
as fit regions which are bounded regularly open sets of finite perimeter and of 
negligible boundary. These postulates enabled the authors to apply a version 
of the Green-Gauss theorem and consider sets that do not necessarily have 
smooth boundaries as bodies in continuum mechanics for which balance laws 
may be written. 

In |30[ 13 lj . Silhavy considered bodies as sets of finite perimeter in a 
bounded open region of M n . The author employs a weak approach in the for- 
mulation of Cauchy's flux theorem. A weakly balanced Cauchy flux of class 
L 1 is shown to be represented by a Borel measurable vector field q of class 
L 1 with a divergence (in the sense of distributions) of class L l . Silhavy's 
approach gives rise to a Borel set Ao of Lebesgue measure zero such that the 
flux vector field q represents the action of the Cauchy flux for any surface 
whose intersections with Ao has Hausdorff area measure zero. The analysis 
presented in Silhavy's work allows for singularities in the flux vector field 
and provides for the first time the concept of almost every surface. In |31| . 
formal definitions of the concepts of almost every body and almost every sur- 
face are given and a weakly balanced Cauchy flux of class LP is represented 
by a measurable vector field of class LP with a divergence of class LP. 

The notions of almost every body and almost every surface are examined 
in [2] and it shown that the Cauchy flux is determined by its action on a 
collection of rectangular planar surfaces with edges parallel to the axes of 
W 1 . A similar extension for Cauchy interaction is presented in |20] , 
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In |28| . a weak formulation of p-grade continuum mechanics, for any in- 
teger p > 1, is presented in the setting of differential manifolds. Configu- 
rations are viewed as C p -embeddings of the body manifold in the physical 
space and forces are viewed as elements of the cotangent bundle to the in- 
finite dimensional configuration manifold of mappings. Forces are shown to 
be represented by measures on the p-th. jet bundle. Such a measure serves 
as a generalization of the p-th order stress. The representation of forces by 
stress measures enables a natural restriction of forces to subbodies. The 
consistency conditions for a such a system of p-th order forces are examined 
in [29]. 

The term fractal was coined in 1975 by Mandelbrot to indicate a highly 
irregular geometric object (see |19j). Mandelbrot's seminal work was the 
beginning of a very large body of research concerning the fractal properties 
of various physical phenomena. A variety of approaches have been suggested 
for the adaptation of fractal objects to branches of mechanics, e.g., \33 \ I34 [ 

ESI ESI Sol SU E5]. 

In [261 127] . Cauchy's flux theory is formulated using Whitney's geomet- 
ric integration theory |39] and new developments by Harrison |12[ [T3| Q31 
I15j . Bodies are viewed as r-dimensional domains of integration in an n- 
dimensional Euclidean space with r < n. A body is identified as an r-chain, 
the limit of a sequence of polyhedral chains with respect to a norm which is 
induced by Cauchy's postulates. Three types of chains are examined: flat, 
sharp and natural chains, such that 

polyhedral C flat chains C sharp chains C natural chains. 

Flat (n — l)-chains may represent the fractal boundaries of bodies and sharp 
chains are shown to represent even less regular (n — l)-dimensional objects. 
Fluxes of a given extensive property are postulated to be (n — l)-cochains, 
i.e., elements of the dual to the Banach space of (n — l)-chains. By the 
duality structure of Whitney's theory, as one allows for less regular domains 
of integration (chains), the resulting fluxes (cochains) become more regular, 
automatically. 

Rough bodies, introduced by Silhavy [32] . are sets whose measure theoretic 
boundaries are fractals in the sense that the outer normal is not defined 
almost everywhere with respect to the (n — l)-Hausdorff measure. 

The present work, describes a framework where the mechanics of bodies 
with fractal boundaries may be studied. Unlike [27], the point of view of 
geometric measure theory as in [4] is mainly adopted. 

The universal body is modeled as an open subset of IR n and bodies are 
modeled as flat chains. In addition to the properties of the class of admissible 
bodies, special attention is given to the study of the kinematics of such bod- 
ies in space. The appropriate class of admissible configurations appears to 
be the set of Lipschitz embeddings. This class enjoys two significant prop- 
erties. Firstly, the set of Lipschitz embeddings of the universal body into 
space is an open subset of the locally convex topological vector space of all 
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Lipschitz mappings of the universal body into space equipped with the Whit- 
ney, or strong, topology. In addition, for Lipschitz mappings there is a well 
defined pushforward action on flat chains, such as those representing bodies. 
Therefore, the images of bodies under the pushforward action induced by a 
Lipschitz embedding preserve their structure and relevant properties (e.g., 
the availability of a generalized Stokes theorem). 

Adopting the point of view that virtual velocities are elements of the tan- 
gent bundle of the configuration manifold, as the configuration space is open 
in the space of Lipschitz mappings, virtual velocities may be identified with 
Lipschitz mappings of the universal body into space. Considering force and 
stress theory, it is noted that forces which are required only to be continuous 
linear functionals relative to the Lipschitz topology, as would be the ana- 
logue of [28], seem to be too irregular for the setting adopted here. In order 
to constitute a consistent force system which is represented by an integrable 
stress fields, balance and weak balance are postulated. It is shown further 
that balance and weak balance are equivalent together to continuity relative 
to the flat norm of chains. 

The paper is constructed as follows. Sections [2H5] contain a short outline 
of the various notions of geometric measure theory which are used in this 
work. Section [2] reviews the notion of differential forms, currents, flat chains 
and cochains. Section [3] presents sets of finite perimeter as well as the corre- 
sponding definitions for bodies and material surfaces as currents. In Section 
U]we discuss some of the properties of locally Lipschitz maps. In particular, 
the image of a flat chain under a Lipschitz mapping is examined. In addi- 
tion, Lipschitz embeddings and the properties of the set they constitute are 
considered. This enables the presentation of a Lipschitz type configuration 
space in Section [6j In Section [5] we discuss the product of locally Lipschitz 
maps and flat chains. This multiplication operation is used in the definition 
of a local virtual velocity. Our main theorem is presented in Section [7] where 
we prove that a system of forces obeying balance and weak balance is equiv- 
alent to a unique n-tuple of flat (n — l)-cochains. Generalized bodies and 
surfaces are introduced in Sections [H Virtual strains, or velocity gradients, 
stresses and a generalized form of the principle of virtual work are presented 
in Sections [9l and PTOl 



2. Review of elements of homological integration theory 

In this section, some of the fundamental concepts form the theory of cur- 
rents in M. n are presented. Throughout, the notation is mainly adopted from 
[U Chapter 4]. The notion of flat forms needed for Wolfe's representation 
theorem, originally presented in Whitney's Geometric Integration Theory 
|39[ Chapter VII], is formulated in this section by the tools of Federer's 
Geometric Measure Theory. 
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Let U be an open set in IR n and V a vector space. The notation D m ([/, V) 
is used for the vector space of smooth, compactly supported ^-valued dif- 
ferential m-forms defined on U and D m (U) is used as an abbreviation for 
T) m (U,M). The notation d<j) is use for the exterior derivative of <j> G T) m {U), 
an element of D m+l {U). The vector space ® m (U) will be endowed with a 
locally convex topology induced by a family of semi-norms [H p. 344] as in 
the theory of distributions. 

A continuous linear functional T : D m (U) — > K is referred to as an Tri- 
dimensional current in U. The collection of all m-dimensional currents de- 
fined on U forms the vector space T) m (U) which is the vector space dual to 
T) m (U). Let T G T> m {U) with m > 1 then dT, the boundary of T is the 
element of D m -i(U) defined by 

dT(4>) =T(d(j)), for all (f) G S) Tn_1 (f7). (2.1) 

The exterior derivative d is a continuous linear map d : T> m (U) — > D m+1 (U). 
Thus, the boundary operation d : £) m +i(£7) -4 D m (U), viewed as the adjoint 
operator to d, is a continuous linear operator on currents. 

As an example of a 0-current in U, let L n , denote the n-dimensional 
Lebesgue measure in M n . Then, the restricted measure L n lU is the 0- 
current defined as 

L n l U ((f>) = [ (f>dL n , for all € ®°(U). (2.2) 
Ju 

Given rj, a Lebesgue integrable m-vector field defined on U, then, L n A rj 
denotes the m-current in U defined by 

L n A 7/(0) = f H V )dL n , for all (/> G D m (U). (2.3) 
./[/ 

The inner product in R n induces an inner product in f\ m M. n and |£| will 
denote the resulting norm of an m-vector £. Given G £) m ([/), for every 
x EU, 4>{x) is an m-covector, and we set 

||</>(x)|| = sup {4>{ x ){Cj I |£l < L ^ is a simple m-vector} . (2.4) 

The comass of </> is defined by 

M(0) = sup \\<t>{x)\\. (2.5) 

For T G £) m (£7) the mass of T is dually defined by 

M (T) = sup {T (<j))\<t>£ T) m (U) , M (4>) < 1} . (2.6) 

An m-dimensional current T is said to be represented by integration if 
there exists a Radon measure /xy and an m-vector valued, /^-measurable 
function, T, with |T(x)| = 1 for //^-almost ah x EU, such that 

T (</>) = / (J)(f)dn T , for all G D m {U). (2.7) 
./[/ 

A sufficient condition for an m-dimensional current, T, to be represented 
by integration is that T is a current of finite mass, i.e., M (T) < oo. An 
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m-current T is said to be locally normal if both T and dT are represented 
by integration and is said to be a normal current if it is locally normal and 
of compact support. The notion of normal currents leads to the definition 

N (T) = M (T) + M (dT) , (2.8) 

and clearly, every T G D m (U) such that N (T) < oo is a normal current. 
The vector space of all m-dimensional normal currents in U is denoted by 
N m (U). For a compact set K of U, set 

N m ,K (U) = N m (U) n {T | spt (T) C K} . (2.9) 

For each compact subset K of U, define Fx, the K-flat semi-norm on 
D m (C/), by 

F K (4>) = sup {||<Kz)||, \\d<f>(x)\\} . (2.10) 

Dually, the K-flat norm for currents T G T) m (U) is given by 

F K (T) = sup {T (0) | F K {(t>) < 1} . (2.11) 

Note that if T G £> m (?7) such that Fr-(T) < oo, then, spt(T) C K. For a 
given compact subset K C U, the set F m ^K(U) is defined as the Fr- closure 
of N m ^K(U) in J) m (U). In addition, set 

^m(CA) =U i? -,^( f/ )' ( 2 - 12 ) 
K 

where the union is taken over all compact subsets K of U. An element in 
F m (U) is referred to as a flat m-chain in U . 

For T G S) m (C7) with spt(T) C K it can be shown that Fk(T) is given by 

Fk{T) = wf{M (T - OS) +M(S)\S G 2) m+ i(*7), spt(S') C . (2.13) 

By taking 5 = we note that 

F K {T)<M{T). (2.14) 

In addition, any element T G F m ^ (U) may be represented by T = R + 95 
where i? G S) m (l7), S G S m+ i(C7), such that spt(i?) C K, spt(5) C i^, and 

F K {T) = M (R) + M {S) . (2.15) 

Flat chains have some desirable properties. We note that the boundary of a 
flat m-chain is a flat (m — l)-chain. Moreover, as Section |4] will show, the 
flat topology is preserved under Lipschitz maps. From a geometric point of 
view the notion of a flat chain may be used to describe objects of irregular 
geometric nature such as the Sierpinski triangle. The following representa- 
tion theorem reveals the measure theoretic regularity characterization of flat 
m- chains. 

Theorem 1. [H Section 4.1.18] Let T be a flat m-chain in U with spt(T) C 
K. Then, for any 5 > and E = {x | dist (K,x) < 5} C U, the current T 
may be represented by 

T = L n An + d (L n A £) , (2.16) 
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such that T] is an L n l U -summable, m-vector field, £ is a L n l V -summable 
(to + l)-vector field and spt {rf) U spt (£) C E. 

A linear functional X denned on F m {U) such that there exists < c < oo 
with X(T) < cFk(T) for any compact K C U and T 6 F m ^K(U), is referred 
to as a /Zai m-cochain. The flat norm of a cochains is given by 



F(X) = sup{X{A) | ^ G Fm(^), ^(i) < 1, K CU}. (2.17) 



By Theorem [H a dual representation for flat cochains is available by flat 
forms which we shall now introduce. 

Given a differentiable mapping u defined on an open set of M n , its deriva- 
tive will be denoted by Du and its partial derivative with respect to the j-ih 
coordinate will be denoted by Dju. For a smooth m-vector field rj in U, the 
divergence divr/ of ij is an (to — l)-vector field in U defined by 



where dxi, i = 1, . . . , n denote the dual base vectors relative to the standard 
basis ej, j = 1, . . . ,n in M. n [H Section 4.1.6]. For an integrable m-form (f> 
in U, the weak exterior derivative of (j) is defined as an (to + 1) form in U 
denoted by dcfi and such that the equality 



holds for all compactly supported, smooth (to + l)-vector fields r\ on U. 
The weak exterior derivative is simply the exterior derivative taken in the 
distributional sense. Note that d(p is uniquely defined up to a set of L n i_ U- 
measure zero, thus, for € D m (U), the relation d(j) = d(p holds L n l [/-almost 
everywhere. 

Differential forms whose components are Lipschitz continuous are referred 
to as sharp m-forms (adopting Whitney's terminology |39^ Section V.10]). 
By Rademacher's theorem, the exterior derivative for sharp m-forms exists 
L n l [/-almost everywhere and the existence of the weak exterior derivative 
follows. Sharp forms are clearly a generalization of the notion of a smooth 
differential form and a further generalization is given by flat forms where the 
Lipschitz continuity is relaxed. 

Definition 2. An m-form (j) in U is said to be flat if 

F(4>) = sup | j (V(r?) + #(6) dL n ^ < oo, (2.20) 

where r/ and £ are respectively m and (m+ 1) compactly supported, L n l £7- 
summable vector fields such that 




(2.18) 



3=1 




(2.19) 




(2.21) 



8 



LIOR FALACH AND REUVEN SEGEV 



It is further observed that for <p, a flat m-form in U, 

F(4>) = e SSS up\U(x)\\,\\d<P(x)\\}. (2.22) 

Alternative definitions for flat forms may be found in |39[ Section IX. 7] 
and [T71 Section 5.5]. 

Remark 3. For <f>, a flat m-form in U, and u, a flat r-form in U, (f) A u is 
a flat (m + r)-form in U. One may use the definition of the weak exterior 
derivative to show that 

d{<j) Alo) = d<j) Awl (-l) m Adu. (2.23) 

The representation theorem of flat cochains is traditionally referred to 
as Wolfe's representation theorem, |39^ Chapter IX], [U Section 4.1.19]. It 
states that any flat m-cochain X in U is represented by a flat m-form denoted 
by D x such that 

X(L n A V + d{L n AO) = [ \D x (ri) + dD x (0]dL n , (2.24) 

Ju 1 J 

for any r\ and £, compactly supported, L n l [7-summable m and (m + 1)- 
vector fields, respectively. It is further noted that the flat norm F(X) for 
the cochain X is given by 

F(X) = esssup \\\D x (x)\l\\dD x (x)\\} = F(D X ). (2.25) 

The coboundary of a flat m-cochain X is defined as the flat (m + l)-cochain 
dX such that 

dX(^) = X(8A), for all ^ G F m ([/), (2.26) 

where it is noted that the same notation is used for the coboundary operator 
and the exterior derivative. The coboundary is the adjoint of the boundary 
operator and thus a continuous linear operator taking flat m-chains to flat 
(to + l)-chains. It follows from the representation theorem of flat chains 
that the flat (to + l)-cochain dX is represented by the flat (m + l)-form 
Ddx = dD x . The last equality is used as the definition of the exterior 
derivative of a flat form in |39[ Section IX. 12]. 

Given a flat TO-cochain X in U and a flat r-cochain Y in U, then, X AY is 
an (m + r)-cochain represented by the flat (m + r)-form D XA y = D x A Dy, 
and for a flat (to, + r)-chain T = L n A i] + d (L n A £), X A Y(T) is defined by 
Equation (|2.24p . Moreover, Equation f|2.23|) implies that 

d{X AY) = dX AY + {-l) m X AdY. (2.27) 

For a flat ?n-cochain X and a flat r-chain T, such that to < r, the interior 
product X jT is defined as a flat (r — m)-chain such that 

IjT(w) = (X Alo) (T), for all w G S r ~ m (?7), (2.28) 

where J Aw is the flat r-cochain represented by the flat r-form D x A oj. 
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3. Sets of finite perimeter, bodies and material surfaces 

In this section we lay down the basic assumptions regarding the material 
universe. Sets of finite perimeter, or Caccioppoli sets, will play a central role 
in the proposed framework. We first recall some of the properties of sets 
of finite perimeter. Extended presentations of the subject may be found in 

plana si]. 

Let U be a Borel set in an open subset of M n and B(x,r) be the ball 
centered at x € W 1 with radius r. Define the U -density of the point x by 

rT , , L n (Ur\B(x,r)) , , 

d(x,U) = lim - V ' , (3.1) 

where the limit exists. The measure theoretic boundary, T (U), of the set U 
is defined by 

r(U) = {x | < d(x,U) < 1}. (3.2) 

Definition 4. A Borel set U in M. n is said to be a set of finite perimeter 
if L n (U) < oo and H n ~ x (T (£/)) < oo, where B n ~ x (r {JJ)) is the (n - 1)- 
Hausdorff measure of T (U). 

Definition [JJ is adopted in |43j as the definition for the class of admissible 
bodies. Several equivalent definitions for a set of finite perimeter may be 
found in the literature. In |44[ Section 5.4.1] a set of finite perimeter is viewed 
as a set whose characteristic function is a function of bounded variation. In 
[H Section 4.5] a set of finite perimeter is viewed as a set which induces a 
locally integral current. In this work, Definition UJ is selected for its intuitive 
geometric interpretation. For a set of finite perimeter, the exterior normal 
v{x) to U exists i5f n_1 -almost everywhere in T(U) thus making a generalized 
version of the Gauss-Green theorem applicable. 

At this point we adopt the framework of |30] for the class of admissible 
bodies and material surfaces. Let B be an open set in W 1 . A body in B is 
denoted by V and is postulated to be a set of finite perimeter in B. Stricly 
speaking, a set of finite perimeter is determined up to a set of L n measure 
zero, thus as a point set, it is not uniquely defined. Formally, each set of finite 
perimeter determines an equivalence class of sets. A unique representation 
of a body is given by the identification of the body V with T-p, an n-current 
in B defined as T v = (L n l V) A e\ A • • • A e n . By Equation (l!H]) . 

T v {oj) = [ u(x) (ei A • • • A e n ) dL n x , for all u € ® n {B). (3.3) 
Jv 

Using the terminology of currents represented by integration, /iy = L n l V 
and T-p = e\ A • • • A e n are the Radon measure and unit n-vector associated 
with the current T-p. 

Objects of dimension (n — 1) for which one can compute the flux will be 
referred to as material surfaces. Formally, a material surface is defined as a 
pair S = (S, v) where S is a Borel subset of B such that for some body V we 
have S C T(V) and v is the exterior normal of V such that v(x) = v-p{x) is 
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defined H n ^almost everywhere on S. Let v*(x) be a the covector defined 
by 

v*(x)(u) = v(x) ■ u, for all u € R n , (3.4) 

and set T$ as the (n — l)-vector 

f s (x) = v*(x) jei A • • • A e n . (3.5) 

It is easy to show that Ts(x) is a unit, simple (n — l)-vector iT"" -almost 
everywhere on S. We use Tg to denote the (n — l)-current in B induced by 
the material surface <S, such that nr s = H™ -1 l5 and Tg(x) are the Radom 
measure and (n — l)-vector associated with Tg, and 

T s (u) = j u{x)(f s {x))dH2~ l , for all u € 2> n_1 (£). (3.6) 

The unit (n — l)-vector Tg(x) is viewed as the natural (n — l)-vector tangent 
to the material surface S. By Equation (|2.3p we may write 

r s = (h*- 1 ^ at s . (3.7) 

Consider the material surface dV = (r(7 3 ),^) naturally induced by the 
body V ■ One has, 

T dP (co) = [ u{x){firp(x))dH^- 1 , 
JT{V) 

(cj(x) jei A • • • A e n ) • v v (x)dH™~ 1 , 

(3.8) 

dw(x) (ei A • • • A e n ) (iL™, 

= (do;) , 
= dT v (u). 

where in the third line above Gauss-Green theorem [H Section 4.5.6] was 
used. Thus, it is noted that T&p = dT-p as expected, and the material 
surface S associated with the body V may be written as 

T s = {dT v )^S. (3.9) 

Since a Radon measure is a Borel regular measure, the current dT-p i_ S is 
well defined for any Borel set <S [U p. 356]. 

For each T v , we observe that M (T v ) = L n (V) and M (dT-p) = ff"" 1 (T(V)) 
correspond to the "volume" of the body and "area" of its boundary, respec- 
tively. By Equation one has N (Tp) = L n (V) + H 71 ' 1 (T(V)) < oo, so 
that the current Tp is a normal n-current in B. The open set B is referred 
to as the universal body and we define the class of admissible bodies, f2g, as 
the collection of all bodies in the universal body B, i.e., 

n B = {T v | V C B,T V = L" L Pe N n (B)}. (3.10) 
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The result obtained in |11| implies that in case B is assumed to be a set of 
finite perimeter, f2g would have the structure of a Boolean algebra and would 
form a material universe in the sense of Noll [23] . In Section EJ a generalized 
class of admissible bodies will be defined for which a requirement that B 
is a bounded set will be sufficient in order to construct a Boolean algebra 
structure. 

The collection of all material surfaces in B will be denoted by dQg, so 
that 

dn B = {t s i t s = (dT v ) l5,t p e n B } • (3.11) 

By the definition of T s it follows that M (T s ) = H n ~ l (s) for each T s £ 

dVLft. Thus T$ is a flat (n — l)-chain of finite mass. The material surfaces 
Ts and Ts> are said to be compatible if there exists a body T-p such that 
Ts = (dTp) l<S and Ts> = (dT-p) i_<S'. The material surfaces Ts and Tg> are 
said to be disjoint if clo(<S) n clo(5') = 0. 

4. LlPSCHITZ MAPPINGS AND LlPSCHITZ CHAINS 

Lipschitz mappings will model configurations of bodies in space. In this 
section we review briefly some of their relevant properties. 

A map T : U V from an open set U C R n to an open set V C M m , is 
said to be a (globally) Lipschitz map if there exists a number c < oo such 
that \J~(x) — J~(y)\ < c\x — y\ for all x, y € U. The Lipschitz constant of T 
is defined by 

fl , = sup gMiM, (41) 

x,j/et/ |y — x\ 

The map T :U — s> V is said to be locally Lipschitz if for every x £ U there 
is some neighborhood f7 x - C f/ of x such that the restricted map J~ \u x is a 
Lipschitz map. 

Let J- : U M. m be a locally Lipschitz map defined on the open set 
U C M. n , then for every K, a compact subset of U, the restricted map T \k 
is globally Lipschitz in the sense that £F,if, the ET-Lipschitz constant of the 
map J- \k, given by 

■W^ir — SU P j j ; l 4 - 2 J 

x,y£K \X — y\ 

is finite. 

4.1. Differential topology of Lipschitz maps. The vector space of lo- 
cally Lipschitz mappings from the open set U C M n to the open set V C M m 
is denoted by £ ([/, V). For a compact subset -fT C U, define the semi- norm 

WH^K = max-dlJ" \k\\oo * ( 4 - 3 ) 

on £ (£/, V), where, 

\K Hoc = SUp (4.4) 
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The vector space £(£/, V) is endowed with the Lipschitz strong topology 
(see [7])- It is the analogue of Whitney's topology (strong topology) for 
the space of differentiable mappings between open sets (see |181 p. 35]) and 
is defined as follows. Given T € £(U,V), for some indexing set A, let 
U = {U\} XeA be an open, locally finite cover of U C W l , V = {V\}^ e ^ 
an open cover of V C M m and JC = {K\} X€A a family of compact subsets 
in U such that K\ C U\ and T(U\) C V\ for all A € A. A neighborhood 
-B £ (J 7 , U, V, 6, fC) of J- in the strong topology is defined by U, V, K, as above 
and a family of positive numbers, 5 = {<5a}a6A' as the collection of all g € 
£(U,V) such that g(K\) C V\ and \\F — g\\^ K < 5\, i.e., 

B z (T,U, V, S, K) = [g e £(U, V) \ g (K x ) c F A; || J" - <7|| £ .^ < 5 X } . 

(4.5) 

Definition 5. A map ip : U — > V, with U C W 1 , V C ffi m , open sets 
such that m > n, is said to be a bi-Lipschitz map if there are numbers 
< c < d < oo, such that [HI p. 78] 

c < \^ x ) ~ V{y)\ < dj for aU x, y E U, x ^ y. (4.6) 
k - y| 

Setting L = max { ^ , d} , 

i<M^M< L , foraU X ,y €U ,X*V, (4.7) 

and in such a case ip is said to be L-bi-Lipschitz. 

The map T : U — > V, where U C W l and V C M m are open sets such that 
m > n, is a Lipschitz immersion if for every x € £/ there is a neighborhood 
U x d U o{ x such that J 7 ^ is a bi-Lipschitz map, i.e., there are < c x < 
d x < oo, and 

|y(y) - 

Cx < ; ; <d x , for all y, z £ U x , y ^ z. (4.8) 

|y - A 

Definition 6. A Lipschitz map (p : U — > V is said to be a Lipschitz embed- 
ding if it is a Lipschitz immersion and a homeomorphism of £7 onto <p(U). 

The following theorem pertaining to the set of Lipschitz embeddings is 
given in [7] and its proof is analogous to the case of differentiable mappings 
as in [Eg p. 36-38]. 

Theorem 7. The set £Em(U, V) is open in £(U, V) with respect to the Lip- 
schitz strong topology 

4.2. Maps of currents induced by Lipschitz maps. Since our objective 
is to represent bodies as currents, and in particular, as flat chains, and since 
we wish to represent configurations as Lipschitz mappings, we exhibit in the 
following the basic properties of the images of currents and chains under 
Lipschitz mappings. 
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Let T be a current on U and for open sets U C M. n and V C W m , let 
F : U — > V be a smooth map whose restriction to spt(T) is a proper map. 
For any r-form uj on V, the map J 7 induces a form J 7 * (cj), the pullback of 
uj by F , defined pointwise by 

(f* (lo) (x)\ ( Vl A • • • A v r ) = (ui (F (x))) (DF(v x ) A • • • A DF(v r )) , (4.9) 

for all Uj, . . . v r G M n . It is observed that since F is proper only on spt(T), 
for a form ui with a compact support, spt(J r *(cj)) need not be compact. 
However, for a real valued function £ defined on U which is compactly sup- 
ported and C(x) = 1 for all x in a neig hborhood of spt (T)Dspt(T*{u)), the 
smooth form (F& (uj) is of compact support. Thus, the pushforward F$, (T) 
of T by F may be defined as the current on V given by 

F# (T) (lo) = T ((F# (oj) \ , for all uj G D r (V) , (4.10) 

for any ^ with the properties given above [HJ Section 2.3]. The definition 
of F# (T) (u) is independent of £ and thus will be omitted in the following. 
The pushforward operation satisfies 

8F # (T) = F#(dT), (4.11) 
spt(J- # T) C J-{spt(T)}. (4.12) 

By a direct calculation one obtains that 

M (F # (T)) < (sup \DF(x)\) M (T) . (4.13) 

\x£K J 



r+1 



Applying Equation (12. 8p it follows that 
N (F # (T)) < N(T) sup \( sup \DF(x)\) Ysup \DF(x)\) \, (4.14) 

[ \x£K J \x&K J J 

and by (j2~T5j) . 

F HK} (F# (T)) < F K (T) sup | (sup |D-F(x)|) , ( sup |DJ-(x) 

[ V^e^ / \x&K 

(4.15) 

where J^jX} is the image of the set if under the map F. 

In case J 7 : U — > V is a locally Lipschitz map, the map F# cannot be 
defined as in the case of smooth maps. However, given any compact K C U, 
for T G F rt x (U), one may define the current F# (T) as a weak limit. Let 
{F T }, t G M + , be a family of smooth approximations of F obtained by 
mollifiers [H Section 4.1.2]. (It is observed that flat chains have compact 
supports so that it is not necessary to require that F is proper.) Set 

F#T(oj) = lim F t# T(lu), for all u G ® r (V). 

T— >0 
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(4.18) 



(4.19) 



The sequence {J~ T # (T)} is a Cauchy sequence with respect to the flat norm 
so that the limit is well defined and one may write 

7#(T) = limJv # (T). (4.16) 

T — 

As a result, the locally Lipschitz map J- : U — > V induces a map of flat 
chains 

J# : F r (U) ->• F r (V) . 

Properties (|4.1ip and f)4.12[) hold for the map Ts, induced by a locally 
Lipschitz map T and 

M (J# (T)) < M (T) (^,spt(T)) r • (4.17) 
It follows that for normal currents 

F#(T) e N r . T{K) (V), for all T G N r , K (U), 

N (J# (T)) < iV(T) sup { (^, S p t (T)) r , (^,s P t(T)) r_1 } , 

and for flat chains 

F#(T) £ F rtHK }(V), for all T G F r , K (U), 

F HK} {T# (T)) < Fk{T) sup { (^, spt (T)) r , (^,s P t(T)) r+1 } • 

See [H Section 4.1.14] and [HJ Section 2.3] for an extended treatment. 

In Whitney's theory, the Lipschitz image of a flat chain A is defined as 
follows |39[ Chapter X]. First, for P = spt(^4) consider a full sequence of sim- 
plicial subdivision {Pi} such that -Pj+i is a simplicial refinement of Pj. Next, 
let {J-i} be a sequence of piecewise affine approximations of the Lipschitz 
map J- such that J~i{v) = F{v) for all vertices v in the simplicial complex 
Pj. The chain T# (A) is defined as the limit in the flat norm of 

T(A) = lim Ti(A). (4.20) 

i— >oo 

Although Whitney's definition of J-jt(A) differs from that of Federer, the 
resulting chains are equivalent. 

For a locally Lipschitz map T : U — > V, and a flat m-cochain X in V, 
let J-# (X) be the flat r-cochain in U defined by the relation 

T* {X){T) = X(T # {T)), for all T G F r (U). (4.21) 

The flat r-cochain J^{X) is represented by the flat r-form J 7 * (Dx), the 
pullback of the flat r-form Dx representing X by the map T . Note that it 
follows from Rademacher's theorem, [4j Section 3.1.6], that DJ- exists L n - 
almost everywhere in U. This does not limit the validity of Equation (14. 9p . 
as a flat form is defined only L n -almost everywhere. 

Consider a locally Lipschitz map T : U — > V from an open set U C M. n 
to an open set V C M m . For a flat n-cochain X in V and a current Tg 
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induced by an L n -summable set B in U, one has 

T* (X) (T B ) = [ F*D x dV\ 
Jb 

= [ D x GF(x)) (DT(x)(ei) A • • • A DF{x){e n )) dL n x , 
Jb 

(4.22) 

D x {F(x)) (ei A • • • A en) Mx)dl%, 




In the last equation the area formula for Lipschitz maps [5J Section 2.1.2] was 
applied and Jj?(x) is the Jacobian determinant of T at x. In case T : U — > V 
is injective with U C M n and 1/ C W 1 , we have 

^ # (X) (T B ) = X (T#T B ) = [ D X (y) dL n y = X (7> {B} ) , (4.23) 

Jr{B} 

thus, J~#Tb = Tjr{ B }. In particular, for a body V, we note that 

T # T v = T r{v] . (4.24) 
For the material surface Tgp, Equation (|4.1ip gives 

F#{T 9V ) = T#(dT v ) = 8T # (T V ) = 8T HV} , 
and for a material surface Ts Equation (13.91) implies that 

T # {T s ) = T ns} . (4.25) 

5. The multiplication of sharp functions and flat chains 

A real valued field over a body V will be represented below by the prod- 
uct of the current T-p and a sharp function — a real valued locally Lipschitz 
mapping. (The terminology is due to Whitney |39t Section V.4].) The space 
of sharp functions will be denoted by C S (U). 

A sharp function (j) € C s (U) defines a flat 0-cochain as on U as follows. 
Let £ be an L n i_ [/-measurable function compactly supported in U. Then, 
L n A £ is a 0-current of finite mass in U as defined in Equation (|2.2p . We set 

a^(L n AO = [ <f>(x)^(x))dL n . (5.1) 
Ju 

For a compactly supported L n i_ U measurable 1-vector field rj, L n A rj is a 
1-current of finite mass in U defined in Equation (|2 . 3[) . Using the existence 
of the weak exterior derivative d<j), L n l [/-almost everywhere, we set 

a^(d(L n At 1 ))= [ d<P(ri(x))dL n , (5.2) 
Ju 

and obtain expressions analogous to Wolfe's representation theorem (Equa- 
tion ([2T24"]) ). Let A G F (U) be a flat 0-chain in U. Applying Theorem [Tj .A 
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may be expressed as A = L n A £ + d (L n A r\) with £ and r\ as denned above. 
Set 

a4A) = a li> (L n At + 8(L n At])), (5.3) 

so that ou defines a continuous, linear function of flat 0-chains. Applying 
Equation (|2.25p we obtain 

F(a ) = sup{|<A(x)|,|#(x)|}. (5.4) 

For .A € F r (U) and € £ s (£/), define the multiplication 0„4 by 0.A = 
a^-iA using the interior product as defined in Equation (|2.28p . That is, 

<j)A{uj) = (a^^A) (w) = (o^/\u)(A), for all u € £ r (t/), (5.5) 

where Aw is the flat r-cochain represented by the flat r-form <f)Auj. Note 
that by Equation ()5.5|) 



spt (<j)A) c spt (0) n spt (A) . (5.6) 

For the boundary of (pA we first note that 

d {4>A) (uj) = 4>A (du) = (a^ A duj) A, for all wer 1 ^). (5.7) 
By Equation ([237|) 

d (o^ Aw) = (da^) Aw + d^A dw, (5.8) 

so that 

9 (w) = (d A w) - (dcty) A a/) A 

= (09.4 - dc^ j .4) (a;) . (5 ' 9) 

Hence we can write 

d (<fiA) = <j)dA — da^ j A. (5.10) 

Remark 8. The multiplication of sharp functions and chains was originally 
defined in |391 Section VII. 1] using the notion of continuous chains which are 
r-vector field approximations of r-chains. 

Proposition 9. Given a sharp function (j), for A € N r ^(U) 

N r , K (<PA) < ( sup \<f>(x)\ + r^jA N r>K (A) , (5.11) 
and for A € F r>K {U) (see [39l p. 208]; 

F r ,K (M) < ( sup \4>(x)\ + (r + 1) Si^K ) F r , x (X) • (5.12) 
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Proof. For A £ N r ^(U) we have 

M( M )= sup 

(c^ A w) (.4) | 



sup 



eD'(ir) M M 



= sup 

u&D*(U) M M 

sup^|| (^(g)g;(x))||M(.4) 

_^ slip , r 

wester) ^ ( w ) 

< sup \4>(x)\M (A) , 



where in the third line we used the representation by integration of A and 
in the fourth line the term sup^g^- |^>(x)| was extracted since spt(*4) C K. 
In order to examine the term M (d (eft A)), we first apply Equation (15.10|) 

M (d {4>A)) < M {4>8A) + M (d a<j) j A) . (5.14) 

For the first term on the right-hand side we have, 

M (cf>dA) = sup \ a 4>^ dA )\ < f sup l<Kx)l \ M {A) (515) 
For the second term, 

\jjjdoi4, Au [f A j dfi A \ 

M (da^ jA)= sup 



we-ar-HU) M(u) 
sup 



< m ,„ su Pa . eK \\d<j>(x) A u(x) \\M (A) 



Luev r - 1 {u) M(w) (5.16) 

r \ sn VxeK \d^(x)\M (u) M (A) 



< sup 



I J M(u) 
r ( sup \d(j)(x)\ ) M{A), 



\xdK 

where in the third line we used the fact that for an £-form uj and a /c-form ui' 
M (u A (J) < ( l + k k J M (u) M (u/) , (5.17) 

as is shown in |5|. 
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One concludes that 

N ((f) A) = M (0 A) + M(d (<j)A)) , 

< sup \4>(x)\M(A) + sup \</>(x)\M (dA) + r£ (f>K M (A), 

xeK x€K (5.I8) 

< (sup \(j)(x)\ + rQ^ K ) N (A) . 

\x£K J 

For a flat r-chain A £ F r x(U) we use the representation given in Equation 
(12J5|) by A = R + OS so that F K (A) = Af (R) + Af (S) . 
We first observe that 



M(da (/> -iS)= sup — > 
uewiu) M (w) 

(da^ A (5) 
SU P I7T~\ ■ 

u&0r(U),spt{w)(ZK M l W J 

, Af (S) M (da^ A u) 



M(oj) 

I 

Af (w) 



"t 1 )M(u;)sup|#(x)|, 



(5.19) 



and conclude that 

Af (da j 5) < (r + 1) £^Af (5) . (5.20) 
Estimating Fk{4>A), one has 



F K (ij>A)=FK(<l>R + <l>dS), 

< F K (<f>R) + (<j)dS) , 

< F K {cj>R) + F K (d a4> j S + d(<t>S)) , 

< Fff {cj>R) + Fr- (dt^ j S) + F K (d(<f>S)) , 

< F K {cj>R) + F K (d a4> jS) + F k ((<PS)) , 

< Af (<f>R) + Af {d a<t> j S) + Af (<f»S) , 

< sup \4>{x)\M (R) + (r + l)£^Af (5) + sup \<f>(x)\M (S) , 

x&K x£K 

< { sup |0(x)| + (r + 1)£^) (Af (J?) + Af (5)) , 
Ue-ff J 

= lsup\0(x)\ + (r + l)£^ K )F(A), 

(5.21) 

where in the third line we used Equation (15.10|) . in the sixth line we used 
Equation (12. 14j) . and in the seventh line we used Equation (I5.20j) . □ 
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The vector space of sharp functions defined on U and valued in R m is 
identified as the space of m-tuples of real valued sharp functions defined on 
U i.e. C s (U,R rn ) = [C a (U)] m . For <f> € C s (U,R m ) and A € F r>K (U) the flat 
r-chain (j)A is viewed as an element of the vector space of (F r ^K{U)) m ', i.e., 
an m-tuple of flat r-chains in U with (4>A) i = 4>iA. 

6. Configuration space and virtual velocities 

Traditionally, a configuration of a body V is viewed as a mapping V —> M. n 
which preserves the basic properties assigned to bodies and material surfaces. 
Guided by our initial definition of a body Tp as a current induced by V , a set 
of finite perimeter in the open set B, a configuration of the body V is defined 
as a mapping K-p £ £Em(^') R n ). To distinguish it from a configuration of the 
universal body to be considered below, such an element, Kp, will be referred 
to as a local configuration. The choice of Lipschitz type configurations is a 
generalization of the traditional choice of C 1 -embeddings usually taken in 
continuum mechanics. 

It is natural therefore to refer to Q-p = £,Em('P,R n ) as the configuration 
space of the body V . Since a body is a compact set, it follows from Theorem[7] 
that Q v is an open subset of the Banach space £(V, R n ) = £ {up {V} , W 1 ). 

For V ,V € the local configurations Kp, K-p/ are said to be compatible 

if 

Kp \~prxp ,= Kp> \~prxp> ■ (6-1) 

Note that the intersection of two sets of finite perimeter is a set of finite 
perimeter, thus, the restricted map may be viewed as the configuration of 
the body Vf\V'. 

A system of compatible configurations AC, IS £1 collection of compatible local 
configurations k = {Kp \ V G ^b}- Clearly, a system of compatible config- 
uration is represented by a unique element of £Em B^)- An element 
k € £ehi (£>, M n ) will be referred to as a global configuration, and the global 
configuration space Q is the collection of all global configurations, i.e., 

Q = £ Em (£,in. (6.2) 

We will view the configuration space as a trivial infinite dimensional differ- 
entiable manifold, specifically, a trivial manifold modeled on a locally convex 
topological vector space as in |21|, Chapter 9]. 

It is noted, in particular, that a Lipschitz embedding is injective and 
the image of a set of a finite perimeter in B is a set of finite perimeter in 
W n . In addition, as Section [4] indicates, Lipschitz mappings are the natural 
morphism in the category of sets of finite perimeters and in the category of 
flat chains. Thus, an element k £ Q preserves the structure of bodies and 
material surfaces as required. That is, every k £ Q induces a map k# of 
flat chains. For any Tp 6 fig, the current k# (T-p) is an element of N n (M n ), 
and for any T$ € d£lg, the current k# (T$) is an (n — l)-chain of finite mass 
in R n . By Equations (j424|) and (j42o]l it follows that k # (T-p) = T k{p} and 
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k# (Tg) = T K is\. Applying Equation fl4. 18j) . one obtains for every T$ £ dflg 
that 

M (k# (T S )) < M (T S ) (£ K)< s) n_1 • (6-3) 
By Equation (|4.17p . for every Tp € fig, 

N(k#(T v )) < N{T V ) sup {(fi^r,^)"- 1 }. (6.4) 

For a global configuration k, let k (fig) denote the collection of images of 
bodies under the configuration k, i.e., 

K(n B ) = {K # (T v )\T P Gn B }. (6.5) 

Similarly, the collection of surfaces at the configuration k is 

k (dn B ) = {^# (t s ) I t s e dn B } • (6.6) 

A global virtual velocity at the configuration n is identified with an ele- 
ment of the tangent space to Q at k. By Theorem [TJ £(B,M. n ) is naturally 
isomorphic to any tangent space to Q. Moreover, k induces an isomorphism 
£(<B, K n ) = £,(K{B},W n ) and an Eulerian virtual velocity is viewed as an 
element of £ (k {B} , R n ). In what follows, we refer to £ (k {B} , ]R n ) as the 
space of global virtual velocities at the configuration k and use the abbre- 
viated notation W K for it. Naturally, an element of W K may be identified 
with an n-tuple of sharp functions defined on n{B}, i.e., using the Whitney 
topology on Z(k{B}), W k = [£(«{£})]" 

Focusing our attention to a particular body V , one may make use of the 
approach of |28| and define a virtual velocity of a body V at a configuration 
K-p E Q-p as an element v-p in the tangent space T Kv Qp. It follows from 
Theorem |7| that one may make the identifications T Kv Qp = £(V,M. n ) = 
£(k p {V},R n ). 

Theorem 10. For every body V , and every Kp £ Q-p, and every k £ Q such 
that k \p= Kp, the restriction mapping 

p P :T K Q^T Kv Qp (6.7) 

is surjective. 

Proof. We recall that Kirszbraun's theorem asserts that a Lipschitz mapping 
/ : A — > IR m defined on a set A C R n may be extended to to a Lipschitz 
function F : M n — > R m having the same Lipschitz constant (see [4j Section 
2.10.43] or pU Section 6.2]). It follows immediately that any v v £ £(V,R n ) 
may be extended to an element v £ £(£>,M n ). □ 

Anticipating the properties of systems of forces to be considered below, 
we wish to provide the collection of restrictions of global virtual velocities to 
the various bodies with a finer structure than that provided by the IHI^^- 
semi-norms. In particular, when considering the restriction v \p of a global 
virtual velocity v to a body V, we wish that the magnitude of the resulting 
object will reflect the mass of V . The local virtual velocity for the body Tp at 
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the configuration k induced by the global virtual velocity v E W K is defined as 
the n-tuple of normal n-currents given by the products vk# (Tp) such that 

[vk # {T v )) i = ViK # (Tp) , for all i = l,...,n. (6.8) 

By Equations f|5. 13[) and (|5.1ip . each component [vk.^ (Tp)]. is a normal 
n-current such that 



M ([vk# (T P )]J < sup \vi(y)\M (k # (Tp)) , 
V&k{P} 

< sup \ Vi (y) (£ K . v ) n \M (Tp) 
y£K{V} 



(6.9) 



and 



N([vk#(Tp)].) < sup \ Vi (y)\ +n£ ViMV} N(k # (T v )) 
V \y^{v} J J 



< 



sup \vi{y)\ +n£ VtMV} 
,y&n{p} I 



(6.10) 



suv{(Z^p) n ,(2 K ,p) n - l }N(T v ) 



In other words, the mapping W K x J7g — s> 33 m (23) given by (v,Tp) i-4 
uft| (Tp) is continuous with respect to both the mass norm and the nor- 
mal norm. 

Similarly, the assignment of a virtual velocity v € W K to a material sur- 
face Tg induces an n-tuple of (n — l)-chains defined by the multiplication 
vk# (T|s). Each component [vk# (Tg)],- is a chain of finite mass and applying 
Equation (|5.13p . one obtains 

M(K # (T 5 )].)< f sup \ Vl (y)\\M(K # (T s )), 
\y&*{s} 



n-l 

< ( sup \ Vi (y)\ | [£ kS ) M(T S ) 
ydK{s} 



(6.11) 

n-l 

.. \ A/r (t„\ 

K,S 



7. Cauchy fluxes 

Alluding to the approach of }28| again, a force on a body V at the config- 
uration Kp € Qp is an element in the dual to the tangent space, T* v Qp. In 
other words, forces on V are elements of the infinite dimensional cotangent 
bundle T*Qp. For gp € T* v Qp, and vp G T Kv Qp, the action gp(vp) is 
interpreted as the virtual power performed by the force gp for the virtual 
velocity vp. It follows immediately that a force on a body V at Kp may 
be identified with a linear continuous functional on the space of Lipschitz 
mappings. Such functionals are quite irregular and will not be considered 
here. 
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Instead, we use in this section the notion of a Cauchy flux at the con- 
figuration real valued function operating on the Cartesian product 
K(dQs) x W K . These impose stricter conditions on the force system and 
resulting stress fields. The conditions to be imposed still imply that for a 
fixed body, a force is a continuous linear functional of the virtual velocities 
of that body. 

A Cauchy flux represents a system of surface forces operating on the ma- 
terial surfaces, or more precisely, their images under k. For a given surface 
and a given virtual velocity field, the value returned by the Cauchy flux 
mapping is interpreted as the virtual power (or virtual work) performed by 
the force acting on the image of the material surface under k for the given 
virtual velocity. 

Definition 11. A Cauchy flux at the configuration k is a mapping of the 
form 

<s> K -.K(dn B )xw K ^R, (7.1) 

such that the following hold. 

Additivity: §> K (-,v) is additive for disjoint compatible material sur- 
faces, i.e., for every k# (Tg) , k# (Tgi) G k(8Qb) compatible and 
disjoint, 

$ K (k# (Tsus>) ,v) = $ k (k# (T s ) ,v) + $ k (/s# (TsO , v) , (7.2) 

holds for every v G W K . 
Linearity: $ K (n# (Ts) , •) is a linear function on W K , i.e., for all a, (3 G 

R and v, v' G W K , 

$ K (k# (T s ) , av + f3v') = a<$> K (k# (T s ) , v) + /3$ K (k# (T s ) , v') (7.3) 

holds for every k# (Ts) G k (8Qb)- 
Let v G W K and k# (T5) G K,(dO,!s), then, by the linearity of the Cauchy 
flux, 

(n \ n 

(t s ) , )= Yl ®- ( K # ( T5 ) ' Vie ^ • ( 7 - 4 ) 
i=l / i=l 

Set («# (T s ),u) = $ K (k # (T 5 ) ,uei) for all it G Z(k{B}), so that is 
naturally viewed as the i-th component of the Cauchy flux at the configura- 
tion k. One has, 

n 

$ K («# (T s ) ,v) = J2$k ( K # ( T s) , Vi) • (7.5) 
i=i 

Balance: There is a number < s < 00 such that for all components 
of the Cauchy flux 

& K (k # (T S ) ,v)<s \\v\\ £tS M (k (T S )) , (7.6) 
for all k # (T s ) G K (8CIb) and v G W R . 
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Weak balance: There is a number < b < oo such that for all com- 
ponents of the Cauchy flux 

K (« # (dT P ) ,v)<b \\v\\ £tV M {k # (Tp)) , (7.7) 

for all k# (Tp) € n (fig) and v £ W K . 

It is observed that from the balance property assumed above, for each 
material surface T$, & K (k# (T$) , •) is continuous. 

Theorem 12. Each component of the Cauchy flux $ K induces a unique flat 
(n — l)-cochain in k{B}. 

Proof. Let o~ n ~ l be an oriented in — l)-simplex in k{B}. Since k {B} is 
open there exists some n-simplex a n in k{B} such that cr n_1 C da n . Since 
k^ 1 {cr n } is a set of finite perimeter in B it follows that a n ~ l £ K (dflg). In 
other words, every oriented (n — l)-simplex in k {B} may be viewed as an 
element of k (dilis). 

In what follows, we use extensions of Lipschitz mappings as implied by 
Kirszbraun's theorem. First, define a real valued function a of (n — 1)- 
simplices. Let u : k {B} — )• K be a locally Lipschitz function in At {0} such 
that u(x) = 1 for x £ <7™ _1 , and we set 

«(a"- 1 )=^(a"- 1 ,u). (7.8) 

The fact that the definition is independent of the choice of u follows from 
condition (|7.6p and will be demonstrated below where a is extended to poly- 
hedral (n — l)-chains. 

Consider a polyhedral (n — l)-chain A = Ylj=i a j°' 1 j ~ 1 * n K {&} sucn that 

•fcr" -1 ! are pairwise disjoint. Define the function u : u/iCT™" 1 — > M by 

u(x) = aj ifarGo-? -1 . (7.9) 

We now apply Kirszbraun's theorem and obtain u : K {B} — > R, a Lipschitz 
extension to u defined on k{B}. By the properties postulated for Cauchy 
fluxes 

%{y j j=1 a--\u) = E*«(°r 1 » fi )=E a i a ( o r 1 )- (7 - 10) 

3=1 3=1 

The function a is now extended to polyhedral (n — l)-chains in k {£>} by 
linearity, i.e., 

a (A) = a[Y,^T l \ =E a i°( T 1 )- (7 ' n) 

Thus, a is a linear functional of polyhedral (n — l)-chains. The value of 
a (A) is independent of any particular extension of it, for given u' , u any two 
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Lipschitz extensions of u, 



(VJ J j=1 a]-\u-u' 



< S ||n-n'L Mj „ n -i M(T n j „ n -i), ( 7 - 12 ) 
= 0. 



From Equation (j7.6[) it follows that 

\a (a^ 1 ) \ < sM (a™" 1 ) , for all a n ~ x G as {B} , (7.13) 
and by Equation (17. 7h . 

\a (da n ) | < bM (a n ) , for all a n G k {£} . (7.14) 
The flat norm of a the functional a is defined by 

F(a) = sup {a(A) \ A is a polyhedarl (n — l)-chain, 

Fa-(^)<1, KCk{£}}, 

and we obtain 



(7.15) 



F(m = max < sup — — sup — — > < max{s,5|. 

We also recall, [2] Section 4.1.23], that polyhedral chains form a dense 
subspace of the space of flat chains, specifically, for every A G F n — i k (M n ), 
a compact subset C C whose interior contain K and e > 0, there is 

and a polyhedral (n — l)-chain A £ supported in C such that 

F C (A-A £ )<e. (7.16) 

Thus, for every flat (re — l)-chain A we have a sequence Aj such that 
lim^j,^ Aj = A. The cochain a is uniquely extended a flat (n — l)-cochain 

F A- 

j— too ^3 



such that for every A = lim^^ A 



V(A) = lim a{Aj). (7.17) 

The foregoing part of the theorem is analogous to |39[ Section V.4]. 

In order to complete the proof we need to show that for Kjt (Ts) G K (dQg) 
and v G C B (k {B}) we obtain * (vk # (T s )) = (k # (T s ) ,v). By d Section 
4.1.17] the class of flat chains of finite mass is the M-closure of normal 
currents. The chain vk# (Ts) is a flat (n — l)-chain of finite mass. Hence, 
the sequence of polyhedral (re— l)-chains {Aj}^ =v converging vk# (Ts) in the 

flat norm, has a convergent subsequence {^4j'}°^_ 1 such that {.A/} converges 
to (Ts) in the flat norm and 

M (vk # (T s )) = limM (A r ) . (7.18) 



ROUGH BODIES 25 



By the definition of a and the balance principle, Equation (17. 6ft the se- 
quence {a (Aj>) } v_ 1 is a Cauchy sequence in R since |a (-4. m ) — a (Afc) j < 
sM(A m -Ak). Hence 

lim a(^,) (T 5 ),i;). (7.19) 

j — >oc 

Since ^ is an extension of a it follows that ^(«4.' ) = ct(A'j) and 
I* (m # (T S )) - K (k# (T 5 ) ,t;)| = I* (m # (T 5 )) - Jim a(^.)l, 



= \V(vk#(T s ))- lim |, 

/-►oo 

= I* (vk# (T S )) - * f lim A A j, 
= \9[vk # (T s ) - lim A,-,) |, 

V j'^oo J 

< max {s, 6} lim i 7 (Is) — .4^) = 



(7.20) 

which completes the proof. □ 



The extension of each flat (re — l)-cochain from C R n to R n is 

done trivially by setting its representing flat (n — l)-form to vanish outside 
k{B}. We conclude that a Cauchy flux <1> K induces a unique n-tuple of flat 
(re — l)-cochains in R n such that 

n 

$ K («# (t 5 ) ( r <5)) > ( 7 - 21 ) 

i=l 

for all u € and k# (T5) E K,(dQig). The inverse implication is provided 

by 

Theorem 13. An n-tuple of flat (n — 1) cochains in R n induces by 

Equation (f7l?7|) unique Cauchy flux 

Proof. For each v £ VF K and k#T$, the Cauchy flux <1> K (k# (T5) , u) will be 
defined by Equation (I7.2ip . and by the components 

K («# (T 5 ) , «<) = ¥ {v iK# {T S )) . (7.22) 

The additivity (IT. 2|) and linearity (IT. 3j) properties clearly hold since *S> 1 is a 
linear function of flat (n— l)-chains. For the Balance (IT. 6[) and weak balance 
(17. 7|) properties, recall that since \t l is a flat (n — l)-cochain, there exists 
C > such that for every flat (re — l)-chain A with support in K, we may 
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write |\& l (.A)| < CFk (A). For the balance property 

|^(k# (T s ),v i )\ = \9 i (v iK# (T S ))\, 

<CF K(s) (v iK #(T s )), 
<CM(v t K#(T s )), [ ■ ) 

<C\\ Vi \\ z§ M{k # (T s )). 

For the weak balance 

|** (« # (dT v ) , Vi ) | = |^ {viK # (OTp)) |, 

< CF K{T) (v iK# (dT v )) , 

= CF K ( V) (d (viK # (Tp)) + dv -jT-p) , 

< C [F Kiv) (d (v iK# (T V ))) + F K[P) (dv _i K# (T P ))] , 

< C [F k(p) (v { k # (T v )) + (dv j k # (Tp))] , 

< C [M (v iK# (Tp)) + M (dv j K# (T P ))] , 



sup \vi(x)\M (k # (Tp)) + n£ VjK (p)M («# (Tp)) 



< C 

<C(n + l)||^|| £ift(p) M(K # (T P )) 



(7.24) 
□ 



Thus, Theorems 1121 and 1 1 3 1 rest at e the point of view presented in |27| that 
the balance and weak-balance assumptions of stress theory may be replaced 
by the requirement that the system of forces is given in terms of an re-tuple 
of flat (n — l)-cochains. 

8. Generalized bodies and Generalized surfaces 

The representation of a Cauchy flux by an n-tuple of flat (n — l)-cochains 
enables the generalization of the class of admissible bodies and the intro- 
duction of a larger class of material surfaces. By a generalized body we 
will mean a subset V of the open set B such that the induced current Tp 
is a flat n-chain in B. Note that the general structure constructed thus far 
holds for generalized bodies. For any configuration k € £Em(£>)^ n ), the 
current k# (Tp) is a flat n-chain in W 1 , and the operations (vk# (dTp)) 
and d^ (vk# (Tp)) are well defined. 

Definition 14. A generalized body is a set V C B such that the induced 
current Tp = L n i_V given by 



Tp(u) = j udL n 
Jp 



is a flat n-chain in B. 
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By [H Section 4.1.24] the current T-p is a rectifiable re-current or an integral 
flat n-chain in B. Moreover, we have 

F (Tp) = M (T f ) = L n (v) . (8.2) 

The above definition of generalized bodies implies that a generalized body 
may be characterized as an n-rectifiable set in B [H Section 3.2.14], or alter- 
natively, as an L n -summable set in B. The class of generalized admissible 
bodies is 

n B = {Tp\V cB,T^eF n (B)}. (8.3) 

As mentioned in Section [3] tig will have the structure of a Boolean algebra if 
B was postulated to be a bounded set. Since N n (B) C F n (B), it is clear that 
C fig. Given T^,,T^, ; £ tig clearly Tj,^, is an element of tig. Contrary 
to the previous definition of bodies, a generalized body needs not be a set of 
finite perimeter. Although V is a bounded set, its measure theoretic bound- 
ary, r('P), may be unbounded in the sense that lf"' _1 (r(7^)) = oo. Generally 
speaking, the boundary of a rectifiable set may not be a rectifiable set. A 
classical example of such a generalized body in R 2 is the Koch snowflake. In 
|32| . such a body is referred to as a rough body. 

Remark 15. It is noted that although every generalized body V induces an 
integral flat n-chain, not every integral flat represents a generalized body. 
However, it seems plausible that a flat n-class, introduced in [42] . is in one 
to one correspondence with the class of generalized bodies. This issue will 
not be considered in this work. 

Considering a generalized surface, we first note that for a generalized body 
Tp, dTj, is a flat (n — l)-chain in B. In addition, the following argument (|6| 
Lemma 2.1]) indicates that the restrictions of flat chains to general Borel 
subsets are not necessarily flat chains. Let H x<s denote the closed half space 
defined by the linear functional A : M. n — > M such that 

H x ,s = {x G R n | \{x) > s} . (8.4) 

Let Tj, € FK, n (B) be a generalized body in B supported in a compact subset 
K of B, so that dT^, is a flat (re — l)-chain, and consider the chain dTp i_ H\ jS . 
One has, 

F K (dT f l H x ,s) = F K (dT f l H X:S + 9 {Tf, l H XjS ) - 8 (T^ l H XjS )) , 

< Fjc {dTj, L H x , s — d (Tp l H XyS )) + F K {8 (T f l H X)S )) , 

< F K (8Tp l ^ AiS d (Frp l H XiS )) + F K {T f l H x>s ) , 

< M (dT^ l H X)S - 8 {T f L H x , s )) + M (T f l H x>s ) . 

(8.5) 
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Since Tp is a chain of finite mass, M (Tp l H\ tS ) < oo. In addition 

/oo 
M (dTp l F A)S - 9 (7> l H x ,s)) ds = M (Tp) , (8.6) 
-oo 

and so we can show that M (dTp l H\ >s — d (Tp l #a,s)) < 00 only for L 1 - 
almost every s£R. 

In order to define a generalized material surface we follow |32| where the 
various properties of flux over fractal boundaries are investigated. 

Definition 16. For a generalized body V, the subset S C T(P) is said to 
be a trace if there exists a set of finite perimeter M such that S = T(V) flM 
and H n - l {T{V) n T (M)) = 0. Each trace S is associated with a unique flat 
(n — l)-chain given by 

T S = dTp nM -dT M ^V. (8.7) 

For each uj S D n ~ 1 (B) we have 

T s (u)= f duj{e\ A • • • A e n )dL n — f uj{f 9M )dH n - 1 , (8.8) 
JvnM JT(M)nf 

where Tqm is defined as in Equation (|3.5p . The set M, of finite perimeter, 
is referred to as the generator of the trace S and it is shown in |32| that S 
depends on M only through the intersection of dTp with M. 
The collection of generalized material surfaces is defined as 

dtl B = {Tj I S is a trace in sj . (8.9) 

We note that by Proposition [9j for all Tg € (9f2# and v € W K , the multiplica- 
tion vk# (Tg) is an n-tuple of flat (n — l)-chains. Thus, by Theorem 1121 the 

Cauchy flux is naturally extended to the Cartesian product W K x k (^diljs^J ■ 

9. VIRTUAL STRAINS AND THE PRINCIPLE OF VIRTUAL WORK 

For Tp £ dtijs and v 6 W K , d (vk# (Tp)) is an n-tuple of flat (n — 1)- 
chains in B. Thus, ^ (d (vk# (Tp))) is a well defined action of an n-tuple of 
flat (n — l)-cochains on an n-tuple of flat (n — 1) chains. Applying Equation 
I5.1UI for each component we obtain 

n n n 

E ( 5 (w ( t p))) = E - E ^ j K # ( t p)) ■ 

i=l i=l i=l 

(9.1) 

Here a Vi is the flat 0-chain defined in Section [5l 

The terms on the right-hand side of the equation above may be inter- 
preted as follows. The term Y^i=i { v i K # ipTp)) * s interpreted as the vir- 
tual power performed by the surface forces for the virtual velocity v on the 
boundary of the body V . Next, for (d (vk # (Tp))) = -d* (vk # (Tp)), 
the n-tuple of flat n-cochains — d^f is viewed as the body force. Thus 
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the term —d$> (vk# \Tp\) is interpreted as the virtual power performed 
by the body forces along the virtual velocity v on the body V . Finally, 
Yli=l {d a vi J K # C^p)) i s interpreted as the virtual power performed by 
the Cauchy flux along the derivative of the virtual velocity v on the body P. 

An internal virtual velocity is viewed as an element upon which the Cauchy 
flux will act. Thus, a generalized internal virtual velocity is defined as an 
n-tuple of flat (n — l)-chains in k{B}. A typical internal virtual velocity 
will be denoted by \ and is viewed as a velocity gradient or a linear strain- 
like entity. Clearly, not every internal virtual velocity is derived from an 
external virtual velocity. Motivated by the above physical interpretation 
and the classical formulation of the principle of virtual work, we introduce 
the kinematic interpolation map 

e : K {tl B ) xW K ^ [F n _! (« (B))] n (9.2) 

such that each component is given by 

(e («# (T f ) ,v)) i = Vi 8k # (T p ) - d (v iK # (Tp)) . (9.3) 

Note that the map e is disjointly additive in the first argument and linear in 
the second. An internal virtual velocity \ ls s &id to be compatible if there 
are V € and v € W K such that 

X = e(K#(T f ) ,v) . (9.4) 

Given a compatible virtual internal velocity % = £ { K # (^p) > v ) we ma Y 
write, 



i=l 

n n 

= J^^i ( V ^# ( 9T v)) ~J2 d ^i ( T r)) ' 

i=l i=l 
n 

= J2da Vi A*i( K# (T^)), (9.5) 

i=l 

and obtain 

* (vk # (dT f )) - cM (vk # (Tp)) = * (e (k # (T f ) ,«)) , (9.6) 

for all € fig and u € W K . We view the last equation as a generalization 
of the principle of virtual power. 



10. Stress 

Applying the representation theorem of flat cochains, a Cauchy flux is 
represented by an n-tuple of flat (n — l)-forms in k{B}. Let denote the 
flat (n — l)-cochain associated with the i-th component of the Cauchy flux. 
Then, D^. will be used to denote its representing flat (n — l)-form. The 
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ra-tuple of flat (n — l)-forms in k {£>} representing the Cauchy flux will be 
denoted by D^> and will be referred to as the Cauchy stress. 

Using the representation theorem for flat forms we obtain an integral 
representation of the principle of virtual power given in Equation (19. 6p . The 
virtual power performed by surface forces is represented by 

n 

n 

(^KAfi))) {T-p) , 

n „ 

V/ d{v i D^ i {K{x))){DK{x){e 1 ) A---ADK(x)(e n ))dL%, (10.1) 
i=i J, P 

n „ 

V / d (v i D^ i (as (x))) (ei A • • • A e n ) J K {x)dL n x . 
~T if 



i=l 



i=l 



Equations (|4.2ip and (|4.9|) were used in the first and second lines. The 
virtual power performed by body forces is represented by 

n 

-J2 d ^ ( T v)) 

i=i 

n 

i=i 

n „ 

= -J2 (vidDvi (k {x)j) (D«(x)(ei) A • • • A D«(x)(e„)) di£, (10.2) 
i=i ^ V 7 

n 

= - V / [VidDiBi (as (x))) (ei A • • • A e n ) J K (x)dL™. 
i=i J -P V 7 

The virtual power performed by internal forces is represented by 

n 

J> # (dc^.A*;) (T^), 

8=1 

^ / (duj A D 9i (as (a;))) (DAs(x)(ei) A • • • A £>Av(x)(e n )) dL™, (10.3) 

8=1 , ' :P 

n „ 

J] / (dvi A D«r ( (as (a;))) (ei A • • • A e n ) J K (x)dL™. 

8=1 ^ 



8=1 

n 
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For k : B — > W 1 , a Lipschitz map, k^^> is an n-tuple of flat (n — 1)- 
cochains in 0. Each cochain K^^i is represented by a flat (n — l)-form 
D K #y. = K^ z D^ i . The associated n-tuple of flat (n — l)-forms, K^Dy is 
identified as the Piola-Kirchhoff stress 

(n*Dy(x)y = J K {x)D^ i (k(x)) . (10.4) 
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